Imaging high-dimensional spatial entanglement with a camera 
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The light produced by parametric down-conversion 
shows strong spatial entanglement that leads to vio- 
lations of EPR criteria for separability. Historically, 
such studies have been performed by scanning a single- 
element, single-photon detector across a detection plane. 
Here we show that modern electron-multiplying CCD 
cameras can measure correlations in both position and 
momentum across a multi-pixel field of view. This capa- 
bility allows us to observe entanglement of around 2500 
spatial states and demonstrate Einstein-Podolsky-Rosen 
(EPR) type correlations by more than two orders of mag- 
nitude. More generally, our work shows that cameras can 
lead to important new capabilities in quantum optics and 
quantum information science. 

Over the last two decades the spatial correlations be- 
tween photon pairs produced by spontaneous paramet- 
ric down conversion (SPDC) have been investigated in 
a number of different configurations, playing an impor- 
tant role in fundamental tests of quantum mechanics and 
in applications for quantum communication and infor- 
mation processing [Ij. The high dimensionality of the 
transverse spatial degrees of freedom (DOF) of the pho- 
ton pairs can be explored either by projections onto a 
discrete basis [2^, such as the Laguerre-Gaussian modes, 
or by using a continuous basis defined by the transverse 
position or momentum of the photons [3]. For instance, 
it has been shown that their transverse position and mo- 
mentum display EPR-type correlations [HIS]. However, 
most measurements to date have relied upon the scanning 
of a single avalanche photodiode or a very small number 
of individual detectors This sequential scanning or 
use of a small number of detectors negates any informa- 
tion capacity advantage in the use of spatial states. In all 
protocols for high-dimensional quantum key distribution 
[H |9], quantum computation [3 and teleportation [10] 
using spatial states, it is essential to perform a full field 
measurement of the photon transverse position, which is 
made possible by using a 2D detector array. 

As with any quantum observable, the transverse posi- 
tion or momentum of single photons can only be mea- 
sured with a limited resolution [11 . In most of the cases, 
spatial entanglement has been detected through measure- 
ments of intensity correlation by scanning the detectors 
in two conjugate planes [H [T2HT5] . By knowing the in- 
tensity correlation distributions, non-separability [121 [14] 
or EPR [H |13l [15] criteria could be tested. In these 



cases, the resolution of the measurements is limited by 
the size of the pinhole or slit used for the scanning. In 
Ref. [16], the authors performed a direct measurement 
of the Schmidt number associated with the transverse 
distribution of the two-photon field, where the intensity 
distribution of only one of the down-converted fields was 
measured in both image and far field planes with the help 
of an intensified CCD. In this case, the resolution of the 
measurement is limited by the size of the camera pixels. 
Sub-shot-noise measurements of correlations in intensity 
fiuctuations in SPDC have already been performed with 
the help of a CCD camera |17]. In another experiment a 
CCD camera was used to measure the complementarity 
between one and two-photon interference from a Young's 
double slit [18 . Nevertheless, noise in CCD cameras pre- 
vents single photon discrimination. 

Recent years have seen a rapid advance in imaging 
technologies. Developments in low- noise electron multi- 
plying CCD (EMCCD) cameras suggest they are capable 
of greater than 90% quantum efficiencies and, if operated 
in a low light environment whilst being maintained at low 
temperatures, can provide single photon sensitivity [T9] . 
EMC CDs have been used in measurements of sub-shot- 
noise correlations of intensity ffuctuations ^Ol [2j as well 
as in intensity correlations [22| in the far-field of the pho- 
ton pairs from SPDC. 

Here we report the observation of spatial correlations 
from SPDC using an EMCCD camera in both the im- 
age plane and far-field of the non-linear crystal showing 
position correlation and momentum anti-correlation, re- 
spectively. In both cases the images are sparse, that is 
the number of recorded photons per pixel is around 0.02. 
We make no post selection of optimum images [23j ; rather 
we obtain correlations by summing over all the recorded 
frames. The use of a 2D array of detectors allows the ex- 
ploration of both spatial DOF of the photons at the same 
time, thereby increasing the number of accessible states 
and the information content per photon [24]. Here we 
use a 201 x 201 array of pixels, to show 2D correlations 
in both position and momentum of the photon pairs. In 
both bases we show correlations of around 50 modes in 
both transverse directions, giving in principle access to 
several thousand entangled spatial states. 
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RESULTS 

Imaging system. The experimental setup used to 
measure spatial correlations in both position and mo- 
mentum is shown in Fig.[l] A 150 mW, high-repetition- 
rate 355 nm laser is attenuated to 2 mW and using a sim- 
ple telescope the Gaussian beam radius is expanded to 
cTp = (0.66 ± 0.05)mm. A 5 mm long /3-Barium Borate 
(BBO) crystal cut for type-I phase matching was angled 
to provide near-collinear output for degenerate down- 
converted photons at 710 nm. The subsequent choice of 
lens configurations was chosen to maintain a constant 
separation (700 mm) between the BBO crystal and the 
EMCCD (Andor iXonS), thus providing an efficient sys- 
tem for manual switching between the image and far-field 
planes. 

As used by other authors [13l [25] , we use a Gaussian 
model to describe the spatial structure of the two photon 
field at the crystal: 
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where N = l/(7rcr_cr+) is a normalization constant, pi = 
(xi^i/i) is the transverse position of photon i {i = 1,2) 
and a± are the standard deviations of the two Gaussians, 
giving the strength of the position and momentum cor- 
relations, cr_ and cr^^, respectively. We define cr+ to be 
equal to the standard deviation of the Gaussian field dis- 
tribution of the pump laser, a^, with wavelength Xp and 
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FIG. 1: Experimental scheme used to measure position 
and momentum correlations, (a) The imaging system 
used to measure position correlations consists of a telescope 
with lenses 100 mm and 250 mm. (b) Momentum correlations 
are obtained with the implementation of three consecutive 
Fourier systems with lenses 50 mm, 100 mm and 200 mm. 



FIG. 2: Position and momentum correlation functions 
in the x' dimension for an increasing number of im- 
ages, (a) Measured position correlations from imaging plane 
scheme, (b) Measured momentum correlations from far-field 
scheme. Blue circles represents the spatial-correlation and 
red squares represent the reference-correlation. The asym- 
metry of the far-field correlation function is the result of a 
non- uniform near collinear degenerate light cone. 



where L is the length of the SPDC crystal and a = 
0.455 is an adjustment constant [26 . Eq. ([T]) is referred 
to as the transverse wave function of the post-selected 
two-photon field of SPDC for a Gaussian pump beam 
[27] , The joint detection probability is 7^(pi,P2) 
|^(Pi5P2)P- According to the Gaussian model and the 
parameters of our system we predict the number of modes 
for joint detections in both position and momentum to 
be (cr+/cr_)^ ^ 3500. 

When measuring position correlations, we use an imag- 
ing system with a magnification of M — 2.5, com- 
prising two lenses of focal lengths 100 mm and 250 mm 
(Fig.jlji). The far-field image is obtained by using a 
composite Fourier system with an effective focal length 
of /e = 100 mm (Fig.[l]3). Given our imaging and far- 
field system we predict a transverse correlation length 
c^pos Mg- ^ 28 /im and amom = fe/{kcFj^) ^ 17 /im, 
where k = kp/ 2 is the wavenumber of the down-converted 
photons. Note that the position correlation length is 
larger than the pixel size of the camera (16 /im) mean- 
ing that the position correlated photon pairs have a high 
probability of being detected in adjacent pixels which 
removes the need to count multiple photons within the 
same pixel. 

A 10 nm bandpass filter centred at 710 nm with a trans- 
mission efficiency r^fiiter = 90% is located immediately 
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in front of the EMCCD. The EMCCD camera used in 
our experiment is a back-ihuminated 512 x 512 array of 
16 X 16 /im^ pixels optimized for visible wavelengths. The 
quantum efficiency of detection is given as the product 
of the efficiency of the camera, filter and lenses. Thus 
providing that either the signal or idler photon is de- 
tected, we predicted a heralding efficiency for detecting 
its entangled partner to he r] ^ 80%. 

Image analysis. A series of 100, 000 images was 
recorded for image and far-field configurations (see Meth- 
ods). The probability distributions V{p-) and V{p-^) 
were calculated respectively by counting the number 
of coincidences as a function of the difference (imag- 
ing plane) or sum (far-field) of the pixel coordinates in 
both x' and y' directions at the plane of the CCD ar- 
ray. Here we define p_ = pi — p2 and p+ = pi + p2 , 
where Pi = {Pxi^PyJ is the transverse momentum of 
photon i. The transformation from coordinates in the 
plane of the CCD array to the crystal is performed us- 
ing the scaling factors 7pos = 1/M for imaging plane 
and 7mom = {kh)/ fe for far-field. Our result is obtained 
by averaging the spatial correlation function over all 
recorded frames. Within each linear dimension this pro- 
cessing reveals correlations between photon pairs, which 
are clearly visible on top of a broader correlation arising 
from uncorrelated events, as shown in Fig. [2] We mea- 
sured the background correlation by performing the same 
analysis but between consecutive frames to give a refer- 
ence correlation. Fig. [2^ shows the measured coincidence 
counts in the image plane as a function of x'_ coordi- 
nates, while coincidence counts measured in the far-field 
as a function of x'_^ are shown in Fig.[2]3. A correlation 
peak in both position and momentum can already be seen 
in a few images and becomes significantly distinguishable 
above noise as more images are summed. 

By subtracting the reference from the spatial corre- 
lation function it is possible to remove unwanted pixel 
defects and obtain a strong correlation in position and 
anti-correlation in momentum as indicated in Fig. [3^ and 
Fig.jSfD respectively. Since we cannot distinguish the de- 
tection of more than one photon on the same pixel, we are 
unable to measure the correlation function for p[ = p'2^ 
for which a value equal to has been assigned. To con- 
firm that these correlations arise between photon pairs 
the neutral density filter used to attenuate the pump 
beam was placed after the BBO crystal to reduce the 
heralding efficiency to approximately 2.0%, whilst en- 
suring the same photon fiux. As shown in Fig.jsj: and 
Fig.|3]i, we do not observe a strong correlation in these 
measurements; however along the y'_ axis in the image 
plane measurements we observe a correlation indicative 
of charge smearing between pixels in the readout direc- 
tion. 

EPR-type correlations. From the thresholded im- 
ages we can calculate the joint probability distributions 
for both the x' and y' coordinates in the image plane 



and far-field, as shown in Fig.|4] Fig. [4^ and Fig.|4]3 show 
the coincidence counts in the planes (xi, X2) and (^1, ^2), 
while Fig.|4]3 and Fig.|4]i show the coincidence counts in 
the planes (j>xx^Vx'2) and {py^^Py^- In the y' direction we 
observe a strong correlation resulting from charge smear- 
ing, thus preventing any meaningful analysis from image 
plane measurements (Fig.[4]3); however the expected anti- 
correlation in the far-field (Fig.[4]i) is still evident. Note 
that background subtraction has been performed. It has 
been shown that EPR-like correlations can be identified 
by violating the inequality [28l [29j 

ALn(^ik2)A^iJp,Jp,J > ^, (3) 

where A^-^(ri|r2) is the minimum inferred variance, de- 
scribing the minimum uncertainty in measuring the vari- 
able T\ conditional on the measurement of variable r2. 
The minimum inferred variance for x-coordinates in the 
far-field measurement is defined as 

ALn(Pxi|P^2) = ^ V{px^)l:^^(j>x^px^)dp^^. (4) 

We define the minimum inferred variance for the x- 
coordinates in the image plane with a similar approach. 
Fitting Gaussians to the joint probabilities in the x- 
direction (Fig. [4^ and Fig.ji]^) we obtain the variances 
shown in Table |lj Substituting these quantities into 
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FIG. 3: Background subtracted auto-correlation func- 
tions in position and momentum for high and low 
heralding efficiency configurations. Auto-correlation 
functions for the image plane (a) and far-field (b) of the 
crystal, summed over 100 000 images and with the measured 
background correlation removed. In (c) and (d) we show the 
same measurements performed with a lower heralding effi- 
ciency. Variances along the x^-direction of the high heralding 
efficiency configuration are indicated. Due to increased coinci- 
dence counts attributed to smearing, the correlation function 
for 2/- = ±1 were set to in (a) and (c) for clarity. 
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FIG. 4: Joint probability distributions for x' and y- 
coordinates. The probability distributions for joint detec- 
tions in both x' and ^^-coordinates in the image plane (a and 
b) and far-field (c and d) are shown respectively Charge 
smearing gives rise to an artificial correlation in the y' direc- 
tion (b and d). 



Eq. ^ we find 

ALJxi|x2)AL,(p,Jp,J = (6.6 ±1.0) X (5) 
A^iJx2|xi)A^iJp,,|p,J = (6.2 ±0.9) X 10-4^2, (6) 

indicating that the transverse DOF of the two photon 
field exhibits EPR non-locality. This result is in good 
agreement with the theoretical prediction (from Eqs. ([T]) 
and Q) of ~ 3 X 10~^h^ and, although performed with 
background subtraction, is an order of magnitude smaller 
than those in Refs. [4, 13^. 

High-dimensional entanglement. From the joint 
distributions in Fig.jlj we can estimate the ratio {cr-^/cr-) 
for both X and y transverse DOF, in the image plane 
(IP) and far-field (FF). We define the maximum num- 
ber of joint detections for our position and momentum 
measurements to be D^^^ = [{o'x'^/<^x'_){<^y'_^/<^y'_)]iP 
and I^mom = [{^x^_/crx^_^){cry^_/cry'_^)]FF' We must also 
account for the proportion of the beam measured by 
the detector array which is calculated by integrating 
the fitted Gaussian functions over the array size. In all 
cases we find this gives access to more than 85% of the 
available states. Due to charge smearing in the image 
plane measurements we must estimate the dimension- 
ality in y based on that of the x dimension, which is 
supported by the circular symmetry of Fig. [3^. Thus, 
the measured dimensionality in position and momen- 
tum is found to be D^os = D^o^ x 0.95^ ^ 3200 and 
^mom = ^mom >< 0-95 X 0.85 ^ 2500, respectively To 
the best of our knowledge this represents the largest di- 
mensionality for any experiment using entangled spatial 
states of photons. 



Variances Experimental Value 



Amin (2:^1^2) 


(119± 15)/im2 


Amin(^2|xi) 


(119± 15)/im2 




(5.5 ±0.4) X IQ-^/i^/im-^ 




(5.2 ±0.4) X IQ-^h^fim-^ 



TABLE I: Variances obtained from Gaussian fittings of back- 
ground subtracted correlation functions shown Fig.|4] 



DISCUSSION 

We have performed a multimode detection of around 
2500 spatially entangled states of photon pairs produced 
by SPDC. We have shown that it is possible to utilize an 
EMCCD camera to measure spatial correlations between 
photon pairs in both the image plane and far-field of 
the down-conversion crystal. After background subtrac- 
tion, we found that the spatial correlations violate EPR 
criterion. We acknowledge that the use of background- 
subtracted measurements do not meet the strict require- 
ments for tests of EPR non-locality. If the background 
was included in the data, it would lead to an overestimate 
of the standard deviations of the correlations. Further- 
more, in order to satisfy non- locality, the measurement 
of each photon must be performed outside the light-cone 
of its entangled partner. This issue could be solved in 
a modified setup in which the two-photons are split and 
sent to two spatially separated cameras. However, our re- 
sults are in close agreement with the theoretical predic- 
tions, providing evidence that EMCCD technology can 
already be used in the characterization of spatial corre- 
lations of photons. 

The ability to access both transverse DOF of photons 
enables a large state space suggesting the use of spatial 
states in quantum communication and quantum infor- 
mation processing. In this case, information should be 
encoded and processed in an alphabet defined as orthog- 
onal spatial states and subsequently read by projecting 
the full transverse field of the single photon on to a com- 
plete set of orthogonal spatial states. Thus, low-noise and 
efficient 2D detector arrays are essential in the use of spa- 
tial states in these applications, for which our work is a 
significant step forward. Future developments in imaging 
technologies will likely result in 2D detector arrays with 
even lower noise and widen their application in quantum 
optics. 

METHODS 

Thresholding. There are several noise sources that affect 
an EMCCD camera. As well as stray photons, the camera is 
subject to dark noise; thermally excited electrons, clock in- 
duced noise resulting from spurious electrons created during 
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charge transfers between pixels, and readout noise generated 
by the readout register. We characterized the level of to- 
tal noise in our camera when operated at the maximum gain 
by obtaining 100, 000 images while the laser is blocked but 
with the shutter on the camera open. Performing an analysis 
similar to that of Ref. [19] we find a Gaussian distribution of 
readout noise, centred at a value 390 e~ with a standard de- 
viation (Jnoise = 6e~, in addition to a long tail that decreases 
exponentially for higher readout values. Our thresholding is 
performed by subtracting the mean background value for each 
pixel and assigning a 1 to pixels with values greater than 
one standard deviation of noise, and a to all pixels hav- 
ing smaller values. This binary thresholding was found to 
reduce the level of noise whilst maximizing the correlation 
signal strength. 

Camera Settings. The camera was operated at — 85°C, 
with a horizontal pixel shift readout rate of 1 MHz, a vertical 
pixel shift every 0.3 /is and a vertical clock amplitude voltage 
of +4 above the default factory setting. The optimum pho- 
ton flux incident upon the camera was chosen to maximize 
the strength of the correlation signal relative to noise sources. 
In our experiment the image-plane measurements were made 
with an exposure time of 0.4 ms, while the far-field measure- 
ments are obtained using an exposure time of 1 ms, in both 
cases giving an average photon flux of approximately 0.02 
photons/pixel, which is equivalent to the level of noise. A 
region of interest measuring 201 x 201 pixels was chosen to 
observe most of the down-converted field, which resulted in a 
frame rate of approximately 5 Hz. 



Electronic address: |matt hew. edg ar @glasgow . ac . uk| 

[1] Walborn, S. P., Monken, C. H., Padua, S., and 
Souto Ribeiro, P. H. Spatial correlations in parametric 
down-conversion. Phys. Rep. 495, 87-139 (2010). 

[2] Dada, A. C, Leach, J., Buller, G. S., Padgett, M. J., 
and Andersson, E. Experimental high-dimensional two- 
photon entanglement and violations of generalized bell 
inequalities. Nat Phys. 7, 677-680 (2011). 

[3] Tasca, D. S., Gomes, R. M., Toscano, P., Souto Ribeiro, 
P. H., and Walborn, S. P. Continuous- variable quantum 
computation with spatial degrees of freedom of photons. 
Phys. Rev. A 83, 052325 (2011). 

[4] Howell, J. C., Bennink, R. S., Bentley, S. J., and Boyd, 
R. W. Realization of the Einstein-Podolsky-Rosen para- 
dox using momentum- and position-entangled photons 
from spontaneous parametric down conversion. Phys. 
Rev. Lett. 92, 210403 (2004). 

[5] D'Angelo, M., Kim, Y.-H., Kulik, S. P., and Shih, Y. 
Identifying entanglement using quantum ghost interfer- 
ence and imaging. Phys. Rev. Lett. 92, 233601 (2004). 

[6] Leach, J., Warburton, R. E., Ireland, D. G., Izdebski, P., 
Barnett, S. M., Yao, A. M., Buller, G. S., and Padgett, 
M. J. Quantum correlations in position, momentum, and 
intermediate bases for a full optical field of view. Phys. 
Rev. A 85, 013827 (2012). 

[7] Almeida, M. P., Walborn, S. P., and Souto Ribeiro, P. H. 



Experimental investigation of quantum key distribution 
with position and momentum of photon pairs. Phys. Rev. 
A 73, 040301(R) (2006). 

[8] Walborn, S. P. and Lemelle, D. S. and Almeida, M. P. 
and Souto Ribeiro, P. H. Quantum Key Distribution 
with Higher-Order Alphabets Using Spatially Encoded 
Qudits. Phys. Rev. Lett.96, 090501 (2006) 

[9] Walborn, S. P. and Lemelle, D. S. and Tasca, D. S. and 
Souto Ribeiro, P. H. Schemes for quantum key distri- 
bution with higher-order alphabets using single-photon 
fractional Fourier optics. Phys. Rev. A77, 062323 (2008) 
[10] Walborn, S. P. and Lemelle, D. S. and Tasca, D. S. and 
Souto Ribeiro, P. H. Quantum teleportation of the an- 
gular spectrum of a single-photon field. Phys. Rev. A76, 
033801 (2007) 

[11] Rudnicki, L., Walborn, S. P., and Toscano, F. Heisen- 
berg uncertainty relation for coarse-grained observables. 
Europhys. Lett. 97, 38003 (2012). 

[12] Tasca, D. S., Walborn, S. P., Souto Ribeiro, P. H., and 
Toscano, F. Detection of transverse entanglement in 
phase space. Phys. Rev. A 78, 010304(R) (2008). 

[13] Tasca, D. S., Walborn, S. P., Souto Ribeiro, P. H., 
Toscano, F., and Pellat-Finet, P. Propagation of trans- 
verse intensity correlations of a two-photon state. Phys. 
Rev. A 79, 033801 (2009). 

[14] Gomes, R. M., Salles, A., Toscano, F., Souto Ribeiro, 
P. H., and Walborn, S. P. Quantum entanglement beyond 
Gaussian criteria. Proc. Nat. Acad. Sci. 106, 21517- 
21520 (2009). 

[15] Walborn, S. P., Salles, A., Gomes, R. M., Toscano, F., 
and Souto Ribeiro, P. H. Revealing hidden einstein- 
podolsky-rosen nonlocality. Phys. Rev. Lett. 106, 130402 
(2011). 

[16] Di Lorenzo Pires, H., Monken, C. H., and van Exter, 
M. P. Direct measurement of transverse- mode entangle- 
ment in two-photon states. Phys. Rev. A 80, 022307 
(2009). 

[17] Jedrkiewicz, O., Jiang, Y.-K., Brambilla, E., Gatti, A., 
Bache, M., Lugiato, L. A., and Di Trapani, P. Detection 
of sub-shot-noise spatial correlation in high-gain para- 
metric down conversion. Phys. Rev. Lett. 93, 243601 
(2004). 

[18] Abouraddy, A. F., Nasr, M. B., Saleh, B. E. A., 
Sergienko, A. V., and Teich, M. C. Demonstration of 
the complementarity of one- and two-photon interference. 
Phys. Rev. A 63, 063803 (2001). 

[19] Lantz, E., Blanchet, J.-L., Furfaro, L., and Devaux, F. 
Multi-imaging and bayesian estimation for photon count- 
ing with emccds. Mon. Not. R. Astron. Soc. 386, 2262- 
2270 (2008). 

[20] Blanchet, J.-L., Devaux, F., Furfaro, L., and Lantz, 
E. Measurement of sub-shot-noise correlations of spa- 
tial fluctuations in the photon-counting regime. Phys. 
Rev. Lett. 101, 233604 (2008). 

[21] Brida, G., Genovese, M., and Ruo Berchera, I. Exper- 
imental realization of sub-shot-noise quantum imaging. 
Nat. Photon. 4, 227-230 (2010). 

[22] Zhang, L., Neves, L., Lundeen, J. S., and Walmsley, I. A. 
A characterization of the single-photon sensitivity of an 
electron multiplying charge- coupled device. J. Phys. B 
42, 114011 (2009). 

[23] Oemrawsingh, S. S. R., van Drunen, W. J., Eliel, E. R., 
and Woerdman, J. P. Two-dimensional wave- vector cor- 
relations in spontaneous parametric downconversion ex- 



6 



plored with an intensified ccd camera. J. Opt. Soc. Am. 
B 19, 102391 (2002). 
[24] Dixon, P. B., Rowland, G. A., Schneeloch, J., and How- 
ell, J. C. Quantum mutual information capacity for 
high-dimensional entangled states. Phys. Rev. Lett. 108, 
143603 (2012). 

[25] Law, C. K. and Eberly, J. H. Analysis and interpreta- 
tion of high transverse entanglement in optical paramet- 
ric down conversion. Phys. Rev. Lett. 92, 127903 (2004). 

[26] Chan, K. W., Torres, J. P., and Eberly, J. H. Transverse 
entanglement migration in hilbert space. Phys. Rev. A 
75, 050101 (2007). 

[27] Monken, C. H., Souto Ribeiro, P. H., and Padua, S. 
Transfer of angular spectrum and image formation in 
spontaneous parametric down-conversion. Phys. Rev. A 
57, 3123-3126 (1998). 

[28] Reid, M. D. Demonstration of the Einstein-Podolsky- 
Rosen paradox using nondegenerate parametric amplifi- 
cation. Phys. Rev. A 40, 913-923 (1989). 

[29] Reid, M. D., Drummond, P. D., Bowen, W. P., Caval- 
canti, E. G., Lam, P. K., Bachor, H. A., Andersen, U. L., 
and Leuchs, G. Colloquium: The Einstein-Podolsky- 



Rosen paradox: From concepts to applications. Rev. 
Mod. Phys. 81, 1727-1751 (2009). 

Acknowledgements 

M.J. P. would like to thank the Royal Society and the 
Wolfson Foundation. M.J.P. and R.W.B. thank the US 
DARPA/DSO InPho program. M.P.E., D.S.T., F.L, 
R.E.W., G.S.B. and M.J.P. acknowledge the financial 
support from the UK EPSRC. J.L., M.A., and R.W.B 
thank Canada's CERC program. The authors would like 
to thank N. Radwell for useful discussions. 

Author contributions 

M.J.P., R.W.B., G.S.B. and J.L. conceived the exper- 
iment. M.P.E., D.S.T., F.L and R.E.W. designed and 
performed the experiment. M.P.E., D.S.T., F.L, J.L., 
M.A. and M.J.P. designed the measurement and analysis 
programs. M.P.E., D.S.T., F.L, J.L. and M.J.P. analysed 
the results. All authors contributed to the writing of the 
manuscript. 



